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Рассмотрено применение стробоскопического отображения Пуанкаре для восстановления невозмущенной 
траектории движения динамической системы в случае возмущения в виде аддитивных помех амплитудного 
характера и при отсутствии априорных сведений о частоте восстанавливаемого сигнала. 
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Введение. Большинство механических систем в своем составе имеют узлы, совершающие враща- 
тельные движения. При стационарном режиме функционирования (неизменные угловые частоты) 
данных узлов такой вид движения характерен тем, что координаты состояния таких систем явля- 
ются периодическими функциями времени или пространственных перемещений. Факторы, вызы- 
вающие периодические вариации координат состояния, разнообразны. При обработке на метал- 
лорежущих станках - это разнообрабатываемость, вызывающая погрешности формы поперечного 
сечения. При рассмотрении качества роторных фрикционных систем - это развитие неравномер- 
ности износа по периоду вращения. К таким же эффектам приводит дисбаланс роторных систем. 
Эта особенность позволяет сделать вывод о том, что проявления возможных дефектов данных 
узлов также будет носить характер, близкий к периодическому. Таким образом, задача диагно- 
стики работоспособности системы может быть сведена к восстановлению периодических состав- 
ляющих в выходном сигнале и анализе причастности тех 
или иных механизмов к их формированию. 
Математическая модель роторной системы. Рассмот- 
рим модель ротора, взаимодействующего со статором че- 
рез трибосреду. Вращение данного ротора осуществляется 
посредством двигателя постоянного тока (рис.1). 

Математическая модель данной системы в форме 
Коши может быть представлена в следующем виде: 
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Рис.1. Модель рассматриваемой системы 
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де С.,С,„,^,,Г, - параметры двигателя постоянного тока; й — коэффициент диссипации; 


®- ускорение свободного падения; ©- угол поворота ротора относительно статора; 
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х»,х, - координаты ротора относительно центра статора; К. (х,х,, а), Р, (хх,,@), 


М (х,х,,@) - силы и момент, формируемые трибосредой при взаимодействии ротора со 


статором (см. ниже). 
Силы и момент РЁ’ (х,х,, 4), Е, (х,х,,%), М (хьх,,%) возникают в результате взаи- 


модействия отдельных элементов поверхности ротора и статора. Причем степень взаимодействия 
зависит от взаимной удаленности этих поверхностей и скорости их движения друг относительно 
друга. Удельная сила взаимодействия может быть разложена на две составляющие: тангенциаль- 
ную и нормальную. Функция, характеризующая зависимость удельной силы взаимодействия от 
взаимной удаленности, называется функцией сближения (рис.2). 


'Для танге 


. 
: 





Рис.2. Пример функций сближения 


Если аппроксимировать функцию сближения полиномом третьего порядка, то зависимость 
тангенциальной и нормальной составляющей удельных сил от взаимной удаленности А можно 
будет записать в следующем виде: 


3 — 
Р(А)=> ^"А;  Р(А)=У МА. 
п=0 г=0 
После интегрирования по взаимодействующей поверхности могут быть получены значе- 
ния: Ё. (х,х,,%); Е, (хьх,,%), М (х,х,, а). Для этого необходимо задать функцию Д(\у). 
Дефекты статора и ротора могут быть представлены как периодические функции Л (у) и 


в. (у) соответственно (рис.3,а). Сдвиг ротора относительно центра обозначим как и. Совокуп- 


ность Л (\), № (\) и и определяют функцию А(\). 
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Рис. 3. Ротор и статор: а — общий случай; 6 -— пример ротора и статора с дефектом в форме импульса 
Функции Л (\) и Л, (\) могут быть представлены в виде ряда Фурье в тригонометриче- 


ской форме с коэффициентами разложения С" и С', соответственно. 
Тогда с учетом вышесказанного силы будут: 
Е. (х,х,, 4) =Р. с0$ф-— РА, 51 ф, 
р, о ых 3шф-+ А, с05ф, 
где ф определяется по формуле: 
м 
ф = агат —, 
Х> 
а силы имеют вид: 


1-(1+®) —1- 


В =т (мс +мс.)+ >С ыс; ое ИИ ) | 


Е К=—ё® 


вет быа кс) 5.6 [рыбы сын СЦ би Кбнны 


где и, - радиус ротора. 
Причем коэффициенты, входящие в формулы сил, могут быть найдены следующим 
образом: 
= +, 
№ = - Л, 


С, е(а-$) +С у ег, К = +]; 
С. = С С", е 9 т КЁ 


С учетом выше определенных коэффициентов момент сопротивления будет: 


М (хьх,, а) = 2 ть | № +№МС, + у С. [ис + у САЗС м 
1=—00 К=Ь—оо 
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Пример функционирования математической модели. Рассмотрим динамику системы для 
случая, когда дефект ротора и статора имеет форму импульса (рис.3, 6). Результаты моделирова- 
ния с учетом шумов, возникающих при измерении, приведены на рис.4 аи 4,6. 

















Рис.4. Результаты моделирования: а — 1- колебания Х,, 2- колебания Х., 


3- колебания угловой частоты вращения; б-— движение ротора в плоскости Х,Х, 


Из рис.4,6 видно, что анализ зашумленных результатов моделирования с целью диагно- 
стирования системы затруднителен. Поэтому необходима предварительная обработка данных. 
Стробоскопическое отображение Пуанкаре. Стробоскопическое отображение Пуанкаре [1] 
для динамических систем, имеющих периодическое решение, определяется как: 

бы =Р(х,,Г), 
где Г - период. 
Иначе говоря, формируется секущая плоскость в фазо- 
вом пространстве (рис. 5). Точки пересечения траекто- 
рии движения с данной плоскостью (с некоторой ее ло- 
кальной областью) образуют множество 


= {Жьх,,*-35х, | ‚ Очевидно, что периодическое дви- 


жение образует множество из одинаковых точек. 
Воздействие внешних помех приводит к варьи- 
рованию точек множества в окрестности точки невоз- 
мущенного движения. Для широкого класса случайных 
возмущений справедливо следующее утверждение: 





Рис.5. Стробоскопическое 1х 
отображение Пуанкаре Хх, = 9%: й 
1=1 


где х, - центральная точка этого множества, причем с ростом и центральная точка стремится к 


точке невозмущенного движения. 
Таким образом, операция восстановления траектории невозмущенного движения может 
быть сведена к нахождению центральных точек при различных смещениях секущей плоскости 
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вдоль траектории движения динамической системы. Назовем это операцией усреднения. Опреде- 
лим операцию усреднения функции АР |. с периодом т следующим образом: 


— )= ву ИИ 


п и+] 
Для определения периода усреднения т воспользуемся теоремой [2]. 


Е(Е+Кт) 


Теорема. Среди множества функций ГР Ор почти наверное, наибольшими и равными ме- 


жду собой значениями средней за период Т’ мощности 
т 
ел, \2 
==|2 (0) 
Го 
обладают те и только те функции, для которых выполняется следующее равенство: т = Ти, где 


иЕМ. 


Так как Г неизвестно, можно использовать Т'>> Г. На рис.6 представлена зависимость 
средней мощности от периода усреднения. 
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Рис.6. Зависимость средней мощности от периода усреднения 
Применим операцию усреднения к ранее полученным результатам моделирования. На 
рис.7,аи 6 показаны восстановленные таким образом функции. 
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а) 6) 
Рис. 7. Результаты усреднения: а- 1- колебания х, 2- колебания х, ‚ 3- колебания угловой частоты вращения; 


6 — усредненная траектория движения ротора в плоскости х, х, 
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Выводы. Построенная модель роторной системы позволяет анализировать влияние пространст- 
венных неоднородностей узлов сопряжения на траекторию движения системы в фазовом про- 
странстве, что в свою очередь, в также оценить отображение пространственной анизотропии на 
наблюдаемые (посредством датчиков виброакустической эмиссии или датчиков угла поворота) 
параметры движения, такие как, например, угловую частоту. Кроме того, предложенный алгоритм 
восстановления, основанный на стробоскопическом отображении Пуанкаре, позволяет выделять 
периодические составляющие в наблюдаемых параметрах, несущие информацию о дефектах уз- 
лов сопряжения. Это имеет важные значения: для обеспечения качества обработки деталей на 
металлорежущих станках и для виброакустической диагностики развивающихся по пространству 
дефектов в роторных системах (например турбинах). 
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М.М. СНУУЕУКО 


АРРЕТСАТТОМ ОЕ ТВОВО$СОРТС РО1МСАВЕ МАР ЕОК КОТОК 
5УЗТЕМ ТМТЕКРАСЕ МООЕ ОЕРЕСТ$ ОТАСМО$1$ 


Тре арр/сайоп о! згоБо$сор/с Ротсаге тар Гог те гесоуегу о! погта!/ ра? о! дупат/с зу&йет тойоп т сазе о! 
реак аа@#ме а!5{игБапсе апа т аб5епсе оРа рпо! даа оп тедиепсу о ге$ога Ее э!дпа! [5 соп$!Еегед. 
Кеу игогаб: 5!дпа! ге${юогайоп, ауегаде Гог {Ве регоа сарасйу, Ротсаге $ гобо$сор/с тар. 
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